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ABSTRACT 


«  The  kill  probability  resulting  from  the  delivery  of  a  salvo  of 

weapons  is  not  a  straightforward  calculation  since  the  aiming  error 
is  likely  to  be  common  to  all  rounds  in  the  salvo.  In  the  absence 
^  of  a  computer  program  or  a  set  of  tables  the  analyst  may  have  to  re¬ 

sort  to  the  binomial  law  l-(l-p)N9  where  p,  the  "single  shot  kill 
probability"  is  computed  on  the  assumption  that  both  aiming  error 
and  round-to-round  error  vary  with  each  round  in  the  salvo.  Use  of 
the  binomial  law  in  the  salvo  case  can  introduce  serious  error. 

As  an  aid  in  determining  weapon  requirements,  comparing  weapon 
systems  effectiveness,  etc,,  salvo  kill  probabilities  against  circular 
targets  are  tabulated  for  a  wide  variety  of  parametric  values  under 
the  following  assumptions;  (1)  one  aims  at  the  center  of  a  target  of 
radius  a  and  fires  a  salvo  of  size  N;  (2)  the  error  in  the  mean  impact 
point  of  the  salvo  from  the  target  center  is  governed  by  a  circular 
normal  density  with  variance  (3)  the  mean  impact  point  is  common 

to  all  rounds  in  the  salvo  but  varies  from  salvo  to  salvo;  (4)  the 
errors  in  the  individual  impact  point  of  shots  within  a  salvo  from 
the  mean  impact  point  are  independently  governed  by  a  circular  normal 
density  with  variance  (5)  the  two  errors  above,  referred  to  as 

the  aiming  error  and  the  round-to-round  error,  respectively,  are 
independent.  The  salvo  kill  probability,  i.e.,  the  probability  that 
at  least  one  round  in  the  salvo  falls  within  the  target,  is  computed 
as  a  function  of  R  =  a/a2,  T  =  0\/02>  and  N  for  R  “  .  1(.  1)3. 0(.  2)5.0, 

T  =  ,1(. 1)3„0(, 2)5.0  and  N  =  1(1)14(2)20,  Various  examples  pertain¬ 
ing  to  the  use  of  the  tables  are  given. 
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I.  INTRODUCTION 


Suppose  one  has  a  circular  target  of  radius  a  with  center  at  the 
origin  of  the  Cartesian  coordinate  system.  Suppose  further  that  a 
salvo  of  N  shots  is  fired  at  the  target  and  that  the  mean  impact  point 
of  the  salvo  is  a  random  point  from  a  circular  normal  distribution  with 
mean  at  the  origin  and  variance  G-^.  Also  assume  that  the  impact  points 
of  the  N  shots  or  rounds  in  the  salvo  are  independent  points  from  a 
circular  normal  distribution  with  mean  equal  to  the  mean  impact  point 


of  the  salvo  and  variance  Gr 


The  problem  is  that  of  finding  the 


probability  that  at  least  one  shot  or  round  in  the  salvo  hits  the 
targets  i.e.,  lies  within  the  target  radius.  Because  of  the  circular 
symmetry  with  respect  to  the  origin,  the  problem  has  recently  been 
coined  the  "axisymmetric"  problem  by  Jarnagin,  reference  Cl],  and  thus 
the.  computed  probabilities  are  referred  to  as  axisymmetric  kill  pro¬ 
babilities  . 


A  typical  salvo  configuration  is  shown  in  Figure  1  for  N  =  6 
where  the  circular  dot  represents  the  mean  impact  point  (drawn  from 
a  circular  normal  distribution  with  mean  vector  ()  and  variance 
and  the  crosses  represent  the  individual  round  impact  points  (drawn 
independently  from  a  circular  normal  distribution  with  mean  equal  to 
the  mean  impact  point  and  variance  G^),  The  key  feature  in  this 
formulation  is  that  it  allows  for  error  in  the  mean  impact  point  of 
the  salvo  (governed  by  ar^)  as  well  as  round-to-round  error. abbht  the 
mean"'. (governed  by  •  The  former  error  is  sometimes  referred  to  as 
the  "aiming"  error  and  is  common  to  all  rounds  in  the  salvo.  For  a 
particular  salvo,  it  results  in  a  definite  bias  in  some  direction  and 
of  some,  magnitude,  but  the  direction  and  magnitude  of  this  error  vary 
from  salvo  to  salvo.  For  example,  if  a  salvo  is  fired  at  map  coordi¬ 
nates  from  a  naval  vessel,  this  error  could  be  due  to  navigational 
and/or  geodetic  inaccuracies.  The  latter  error,  assumed  independent 
of  the  first,  is  sometimes  referred  to  as  the  "round-to-round"  or 
"random"  error  as  it  varies  from  round  to  round  within  the  salvo. 

For  example,  it  could  be  due  to  ballistic  characteristics  peculiar  to 
individual  rounds. 

The  solution  to  the  axisymmetric  problem  has  been  given  by 
Jarnagin,  reference  [ l] ,  and  the  purpose  of  this  report  is  (1)  to 
briefly  discuss  Jarnagin's  solution  with  respect  to  normalization 
of  the  parameters  involved  to  nondimens ional  units  and  (2)  to  set 
forth  an  extensive  tabulation  of  axisymmetric  kill  probabilities. 
While,  to  the  best  of  this  author's  knowledge,  this  is  the  first  table 
of  this  kind  to  be  published  for  circular  targets,  a  table  of  salvo 


1 


kill  probabilities  for  square  targets  was  published  in  1954  by  the 
National  Bureau  of  Standards,  reference  [ 2 J .  The  axisymmetric  case 
under  discussion  is  more  realistic  for  a  nuclear  attack  against  hard 
sites  where  a  "circular  target"  is  synonymous  with  an  assumed  "circular 
weapons  effectiveness,"  and  the  target  is  simply  a  point. 


II.  DISCUSSION 


Jarnagin's  formulation  of  the  expression  for  the  axisymmetric  kill 
probability  is  given  in  detail  in  reference  Cl].  However,  due  to  its 
brevity  and  for  the  sake  of  completeness  in  this  report,  it  will  be 
repeated  below.  Letting  the  circular  target  have  a  radius  a  and  be 
centered  at  the  Cartesian  origin,  consider  a  small  element  of  area  dxdy 
in  the  real  plane.  Under  the  circular  normal  assumption  given  in  the 
introduction,  the  probability  that  the  mean  impact  point  of  the  salvo 
is  somewhere  in  this  element  of  area  is 


2 
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(1) 


2 

2naj_ 


-(x2+y2)/(2a12) 
e  dxdy 


or  in  polar  coordinates  is 


-r2/(2CTl2) 

r  drd9 


(2) 


Given  that,  the  mean  impact  point  of  the  salvo  is  in  the  given  element 
dxdy,  the  probability  that  a  particular  shot  hits  the  target  is  simply 
P (a/cT2 jr/c^)  where  P(R,d)  represents  the  well  known  circular  coverage 
function  which  is  the  integral  of  a  circular  normal  distribution  of 
unit  standard  deviation  over  a  circle  of  radius  R  offset  from  the 
origin  by  a  distance  d.  A  discussion  of  this  function  and  an  efficient 
program  for  its  computation  are  given  by  DiDonato  and  Jarnagin  in 
reference  [3].  The  conditional  probability  (conditional  on  the  mean 
impact  point)  that  a  particular  shot  misses  the  target  is  then 
l-P(a/o‘23r/a2) ,  and  the  conditional  probability  that  all  N  shots  in 
the  salvo  miss  the  target  is 


[l-P(a/a2,r/a2)]N 


(3) 


due  to  the  assumed  independence  of  the  shots  within  the  salvo.  Hence, 
the  probability  that  the  mean  impact  point  of  the  salvo  is  in  the  element 
of  area  dxdy  and  that  all  N  shots  in  the  salvo  miss  the  target  is 


m  r  -r2/ (2ai 2) 
[l-P(a/CT2,r/a2)]  — — —  e 

2nai2 


drd9. 


(4) 


Since  the  mean  impact  point  of  the  salvo  can  be  anywhere  in  the  planes, 
the  probability  that  all  N  shots  miss  the  target,  regardless  of  the 
position  of  the  mean  impact  point  of  the  salvo,  is  given  by  integrating 
this  expression  over  the  entire  plane,  i.e.. 


3 


2tt 


CO 


r  drdQ  . 


(5) 


/  / 
9=0  r=0 


[l-P(a/a25r/a2)] 


N 


-r2/(2CTl2) 


The  probability  A  that  the  target  is  hit  by  at  least  one  of  the  N  shots 
in  the  salvo  is  clearly  one  minus  the  expression  in  (5)  or 


A  =  1 


2tt  00 


[l-P(a/CT2,r/a2)] 


N 


2  ? 

-rH2oi  ) 

e  r  drd9  . 


(6) 


t 


Since  the  integrand  in  (6)  is  independent  of  0,  one  can  integrate  with 
respect  to  0  and  obtain  the  result 


A  =  1 


[l-P(a/s2,r/a2)] 


N 


2  2 

-r  /<2af) 

e 


r  dr  . 


(7) 


As  aforementioned,  the  above  formulation  is  given  in  reference  [1], 

An  alternative  formulation,  leading  to  the  same  solution,  can  be 
obtained  by  writing  the  conditional  probability  of  at  least  one  hit  in 
density  form  and  using  the  classical  probability  rules  pertaining  to 
conditional,  joint,  and  marginal  probability  densities.  This  is  accomp¬ 
lished  by  defining  the  discrete  random  variable  z  to  be  the  number  of 
times  the  target  receives  at  least  one  hit  in  a  salvo  of  size  N„  Clearly, 
z  can  take  on  only  the  values  0  or  1,  and  the  conditional  density  of  z, 
given  a  mean  impact  point  (r,0),  can  be  expressed  as 


f(z|r,0)  =  {l-[l“P(a/a2,r/CT2)]N}Z{[l-P(a/a2,r/cT2)]N}1  Z  (8) 
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for  z  =  0,  1  ;  0  <  r  <  “  ;  0<9<2rr„  One  can  now  obtain  the  joint 
density  of  all  three  variables,  say  g(z,r,0),  by  multiplying  f(z|r,0) 
by  the  joint  density  of  (r,9),  say  h(r,0)o  By  assumption, 

» 

r  -r2/ (2a,2) 

t  h(r,9)  =  - - g  e  <9) 

2TTa! 


so  that 


g(zs,rs9)  =  f (z|  r,9)h(r,0) 


(10) 


for  z  =  0,l;0<r<o:>  ;0<9<2rr  where  f(z|r,9)  and  h(r,9)  are 

given  by  (8)  and  (9),  respectively.  The  marginal  density  of  z,  say 
p(z),  can  now  be  found  by  integrating  g(z,r,9)  over  the  entire  range 


of  both  r  and  9,  that  is 

2tt  » 

P(z)  =  j  J 
0  0 

However,  since  only  the  probability 
interest,  it  is  found  to  be 


g(z,r,9)drd0  .  (11) 


f  at  least  one  hit,  p(l),  is  of 


Pd) 


2rr 


0 


I 


{l-[l-P(a/o,2,r/a2)]N} 


2  ,  2 
r  -r  /(2a1  ) 


drd9 


t 

j 


OO 


0 


[l-P(a/a2,r/a2)] 


N 


r 

2ttct12 


2  2 
-r  /  (2cj|  ) 

e  drd9 
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If  N 

=  1  -  2  [l-P(a/CT2,r/ff2>]  e 

CTl  J 


=r2/( 2a!2) 


r  dr 


(12) 


« 


1 

which  is  precisely  the  form  of  A  in  equation  (7). 

Examination  of  equation  (7)  (or  equivalently  (12))  reveals  that  A 
is  a  function  of  four  paramet:erss  namely  a,  C29  and  N„  Hences  A  can 
be  expressed  in  functional  form  as  A(asC7pCF2sN) .  As  a  function  of  four 
parameters ,  any  tabulation  of  A  poses  difficulties  since  to  cover  a 
practical  range  for  all  parameters  would  require  a  very  long  table  in¬ 
deed.  Fortunatelyj,  the  problem  can  be  circumvented  by  normalizing  the 
parameters  to  nond i mens i anal  units  with  the  effect  of  eliminating  one 
of  the  parameters.  It  was  chosen  to  normalize,  with  respect  to  C2s  and 
it  can  be  easily  shown  that 


A(a3alsCT2  9N)  =  A(a/CT2sO'1/a2  s  1SN)  , 


(13) 


that  iss  that  A  is  a  function  of  only  three  parameters:  the  normalized 
target  radius ,  the  normalized  ffp  and  N.  This  is  shown  by  making  the 
transformation  y  =  r/a2  in  (7)  and  letting  R  =  a/a2  and  T  =  0^/02  which 
leads  to 


„  -y2/(2t2> 

[l-P(R,y)]  e  y  dy  (14) 


which  is  a  function  of  only  R,  T,  and  N.  In  the  tabulation  which 
follows 5  A  in  expression  (14)  will  be  designated  as  simply  A(R,TjN). 


III.  COMPUTATION 

Unfortunate ly 3  the  integral  in  expression  (14)  cannot  be  expressed 
in  closed  form  since  P(R,y)  itself  is  a  double  integral  which  cannot  be 
expressed  in  closed  form.  However ,  due  to  an  efficient  algorithm  devised 
by  DiDonato  and  Jarnagin  for  evaluating  P(R,y),  reference  13],  the 
integral  in  (14)  lends  itself  well  to  numerical  integration.  One  must 
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first,  however,  find  a  large  finite  upper  limit  for  the  integral  which 
will  result  in  negligible  error.  This  was  found  by  Jarnagin,  reference 
[  l]  ,  for  the  integral  in  equation  (7),  and  his  same  procedure  can  be 
applied  to  the  integral  in  expression  (14).  Suppose  one  integrates 
from  zero  to  cT  vice  zero  to  co ,  The  resulting  error  E  is 


f 


/ 


2  2 

N  -y  / (2t2) 

[l-P(R,y>]  e 


y  dy 


2  2 
-y  /(2TZ) 

ye  dy 


(15) 


since  [l-P(R,y)]  is  always  between  zero  and  one.  This  is  the  same 
result  obtained  by  Jarnagin  except  his  c  was  in  units  of  c  and  the 
above  c  is  in  units  of  T  =cr]_AT2»  For  all  computations  displayed  in 
this  table,  c  was  taken  as  5  which  provided  an  error  of  less  than 
.00000373. 

Using  an  upper  limit  of  5T  on  the  integral  in  (14),  the  integration 
was  performed  on  the  CDC  6700  digital  computer  at  the  Naval  Weapons 
Laboratory  using  an  iterative  procedure  based  on  the  parabolic  (or 
Simpson's)  rule.  This  requires  no  elaboration  except  to  say  that  the 
interval  length  was  halved  in  successive  iterations,  and  the  integration 
tolerance  was  taken  as  .00001,  that  is,  when  the  maximum  absolute  dif~ 
ference  between  three  successive  approximations  for  the  integral  fell 
within  this  tolerance,  the  iterative  procedure  was  terminated.  This 
tolerance,  and  a  c=5  are  more  than  adequate  to  claim  the  four  decimal 
place  accuracy  to  which  the  table  entries  are  printed. 

The  computer  time  for  the  entire  table  of  27,200  entries  was  3172 
seconds  or  an  average.  117  milliseconds  per  entry. 


» 


’ll 
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IV.  DESCRIPTION  OF  TABLE 


The  table,  presented  in  Appendix  As  provides  values  of  A(R,TSN) 
for  R  ■  .1(. 1)3. 0(. 2)5.0,  T  =  . 1( . 1)3 .0( .2)5.0,  and  N  =  1(1)14(2)20 
where  R  =  a/cr2s>  T  =  a\/c 2,  and  N  =  the  salvo  size.  For  each  value  of  * 

Ns  the  first  page  contains  values  of  A  for  all  40  values  of  R  and 
values  of  T  from  .1  to  1.4,  the  second  contains  values  of  A  for  all 
40  values  of  R  and  values  of  T  from  1.5  to  2.8,  and  the  third  contains  * 

values  of  A  for  all  40  values  of  R  and  the  remaining  values  of  T. 

Hence,  three  pages  are  required  for  each  value  of  N  which  results  in 
a  total  of  51  pages.  Examples  of  the  use  of  the  table  are  given  below. 

EXAMPLE  1 

Suppose  one  has  a  circular  target  of  radius  20  feet,  an  aiming 
standard  deviation  of  100  feet,  and  a  round-to-round  standard  deviation 
of  50  feet.  Suppose  further  that  a  hit  on  the  target  provides  a  kill 
and  that:  a  salvo  of  size  12  is  fired  at  the  target.  What,  is  the  pro¬ 
bability  of  kill?  The.  probability  of  kill  is  simply  the  probability 
of  at  least  one  hit  in  a  salvo  of  size  12  and  can  be  found  by  entering 
the  table  with 


R  =  a/<72  =  20/50  =  .4 
T  =  a1/o2  =  100/50  =  2 
N  =  12 


One  finds  A(.4,2,12)  =  .1526. 

EXAMPLE  2 

Suppose  one  has  a  point  target  with  negligible  radius  and  is 
firing  projectiles  with  an  effectiveness  radius  of  100  feet.  Suppose 
also  that  the  aiming  standard  deviation  is  200  feet,  the  round-to~ 
round  standard  deviation  is  50  feet,  and  the  salvo  size  is  5.  What 
is  the  probability  of  kill,  i.e.,  the  probability  that  at  least  one 
projectile  falls  within  100  feet  of  the  target?  Here  one  uses  the 
weapon's  effectiveness  radius  as  the  target  radius  and  enters  the  table 
with 


R  »  a /a 2  =  100/50  =  2 

* 

T  =  oxh2  =  200/50  =  4 
N  =  5. 


One  finds  A(2,4,5)  =  .2426. 
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EXAMPLE  3 


Suppose  one  has  a  point  target  of  negligible  radius  and  is  attack 
ing  with  nuclear  warheads  with  an  effectiveness  radius  of  500  yards. 
Suppose  the.  aiming  CEP,  CEPa<,  is  500  yards ,  the  round-to-round  CEP, 
CEPr,  is  250  yards,  and  the  number  of  warheads  is  three.  What  is  the 
probability  of  kill?  To  use  the  table  in  this  report  to  answer  the 
above,  one  must  perform  several  preliminary  computations.  First, 
the  CEP*s  must  be  converted  to  standard  deviations  via  the  relation¬ 
ship  cr  =  CEP/1.1774.  Second,  in  nuclear  damage  assessment,  the 
round -to -round  standard  deviation  is  usually  inflated  through  the 
relationship 


cr2  =  (CEPr/1. 1774)2  +  (ka)2 


where  a  is  the  weapon  effective  radius  and  k  is  a  parameter  reflecting 
the  uncertainty  of  the  weapon  effect.  T|he  parameter  k  is  usually 
assigned  a  value  of  .2  or  .3,  reference  [4].  Suppose  it  is  .2  in 
this  example.  Then  one  enters  the  table  with 


R  =  a/0'2  =  a/  >/(CEPr/1. 1774) 2  +  (72a)2 


=  500/  \/(250/l. 1774)2  +  [(.2) (500) 


=  500/234.70  =  2.13 


T  =  0^2  =  (CEPa/1. 1774)  /234. 70 


=  (500/1. 1774) /234. 70  =  1.81 


N  =  3. 
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One  finds  A(2.1, 1.9,3)  =  .6145  and  A(2.2,1.8,3)  =  .6759  so  that  the 
probability  of  kill  lies  somewhere  between  these  values.  Should  one 
require  more  information  than  this,  he  could  find  an  approximation  to 
A(2 . 13, 1.81,3)  by  using  the  tabular  values  for  arguments  around 
R  =  2.13  and  T  =  1.81  in  one  of  the  many  bivariate  interpolation 
formulas . 

EXAMPLE  4 

Suppose  one  has  a  point  target  with  negligible  radius  and  is 
firing  projectiles  with  an  effectiveness  radius  of  220  feet.  If 
Ol  =  210  feet  and  erg  =  100  feet,  how  large  must  the  salvo  be  to  insure 
a  target  kill  at  the  .70  probability  level?  Here  one  scans  the  table 
for  the  correct  value  of  N  using 


R  =  a /a 2  =  220/100  =  2.2 


T  =  a  1/o,2  =  210/100  =  2.1 


One  finds  A(2.2,2.1,6)  =  .6836  and  A(2.2,2.1,7)  =  .7033.  Hence,  a 
salvo  of  size  N=7  would  be  minimum  to  insure  a  target  kill  at  .70 
probability  level. 


IV.  CONCLUSION 


It  is  hoped  that  the  table  contained  herein  will  fill  one  of  the 
many  voids  currently  existing  in  weapons  systems  analysis.  The  pre¬ 
vious  nonexistence  of  tabular  salvo  kill  probabilities  for  circular 
targets  has  resulted  in  the  use  of  some  rather  crude  approximations 
by  weapons  systems  analysts.  For  example,  it  is  not  uncommon  to  see 
salvo  kill  probabilities  computed  (incorrectly)  by  using  a  common 

2  2  1/2 

standard  deviation  c  =  (a-^  +  a2  ) 


in  the  cumulative  radial  distri¬ 


bution  to  obtain 


p  =  1  -  e 


-a2/ (2cr2) 


(16) 
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The  value  of  p  in  (16)  is  then  used  as  a  "single  shot  kill  probability" 
and  applied  to  the  binomial  law  to  obtain  a  salvo  kill  probability 
equal  to 


1  - 


d-p) 


N 


(17) 


This  procedure  can  result  in  serious  error.  For  instances  consider 
EXAMPLE  4  in  the  last  section  where  it  was  determined  that  a  salvo  of 
size  N=7  was  minimum  for  a  target  kill  at  the  .70  probability  level. 
Employing  the  above  procedure,  one  finds  02  =  (210)2  +  (qoO)2  =  54100 
and 


-(220)  2 /[  2(54100)3 

p  =  1  -  e 


=  .3607. 


Setting  the  expression  in  (17)  equal  to  .70  with  p  =  .3607  results  in 
N=2.7  or  a  salvo  of  size  3.'  This  is  a  serious  departure  from  the 
correct  value  of  N=7  -  an  error  of  57%  in  salvo  size. 

The  above  procedure  is  valid  only  for  the  case  where  N=l.  In 
fact,  it  can  be  shown  that  p  in  (16)  is  identically  equal  to  A(R,T,1) 
so  that  for  this  special  case,  one  can  express  A  in  closed  form.  How¬ 
ever,  the  above  procedure  is  not  valid  for  N  >  1,  and  one  must  resort 
to  numerical  integration.  The  reason  it  is  not  valid  is  that  in  using 
p  as  the  so-called  "single  shot  kill  probability"  in  (17)  one  is  assuming 
that  the  aiming  error  as  well  as  the  round-to-round  error  varies  for 
each  round  in  the  salvo.  This  is  contrary  to  the  assumption  that  the 
aiming  error  is  common  to  all  rounds  in  the  salvo  but  varies  from 
salvo-to-salvo  while  the  round-to-round  error  varies  from  round  to  round 
within  the  salvo. 

The  selected  values  of  the  parameter  R,  T,  and  N  for  which  A  is 
tabulated  should  be  sufficient  for  most  practical  applications.  However, 
should  the  reader  have  need  of  A  for  extended  values  of  the  parameters 
or  for  intermediate  values  for  which  interpolation  is  not  adequate  a 
FORTRAN  listing  of  the  AKP  (Axisymmetric  Kill  Probability)  program 
is  shown  in  Appendix  B.  The  version  shown  is  not  normalized  with 
respect  to  s0  that  A  is  computed  as  a  function  of  a and  N 
as  shown  in  equation  (7). 
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APPENDIX  A 


Table  of  A(RST,N) 


fl(S,T,N)  =  ERGE ( AT  LEAST  ONE  HIT)  FCR 
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ter  of  a  target  of  radius  a  and  fires  a  salvo  of  size  N;  (2)  the  error  in  the  mean 
impact  point  of  the  salvo  from  the  target  center  is  governed  by  a  circular  normal  den¬ 
sity  with  variance  o  (3)  the  mean  impact  point  is  common  to  all  rounds  in  the  salvo 
but  varies  from  salvo  to  salvo;  (4)  the  errors  in  the  individaul  impact  point  of  shots 
within  a  salvo  from  the  mean  impact  point  are  independently  governed  by  a  circular  nor¬ 
mal  density  with  variance  a??;  (5)  the  two  errors  above, referred  to  as  the  aiming  error 
and  the  round-to-round  error,  respectively,  are  independent.  The  salvo  kill  probabili¬ 
ty,  i.e.,  the  probability  that  at  least  one  round  in  the  salvo  falls  within  the  target, 
is  computed  as  a  function  of  R  =  a/o;,  T  =  o ,/«;•,  and  N  for  R  =  .1  (.1)3. 0(. 2)5.0, 

T  =  ,1 (  1 ) 3 - 0 ( .2)5.0  and  N  =  1(1)14(2)20.  Various  examples  pertaining  to  the  use  of 
the  tables  are  given. 
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